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Abstract 

We prove a formula conjectured in 14 for the free energy density 
of a directed polymer in a Brownian environment in 1 + 1 dimensions. 
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1 Introduction 

Let B^^\B^'^\ ... be independent standard one-dimensional Brownian mo- 



tions. Denote the increments of B^^^ by B^^^,. = B^^^ — B^^^ 



For /? e M set 

ZniP) = [ ds,...dsn-iexp{P{B[^^^^^ + ... + B^^l^^^^)}. (1) 

This is the partition function for a continuous model of a directed polymer 
in a Brownian environment in 1 + 1 dimensions. In the paper ^Ij, using 
queueing-theoretic ideas in the context of geometric functionals of Brown- 
ian motion, certain limiting results were obtained which led the authors to 
conjecture an explicit formula for the free energy density 

lim -logZ„(/3), (2) 

ra— »oo n 

namely that it should be given by, almost surely, 

/(„) = {-(-*>•(-« -2'°^ I"! :;^f» (3) 

where ^(m) = T'{m)/T{m) is the restriction of the digamma function to 
(0, oo), and (—'I')* is the convex dual of the function The aim of this 
paper is to give a rigorous proof of this conjecture. The proof uses tools 
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from large deviation theory. As a corollary we give a new proof that c > 2, 
where 



lim —Ln{n) = c a.s. 

n— >oo n 



(4) 



and 



Ln{t) 



sup 

0<si<...<s„_i<t 



^(1) 
(0,si) 



+ ... + B, 



(Sn-lv 



It was proved using direct methods in that c = 2, where the authors 
also describe how the result may be deduced from the theory of random 
matrices. 

The directed polymer model we have discussed here, and for which we 
have computed the free energy density, is a continuous version of the classical 
two-dimensional directed polymer, where it is not known how to compute 
the free energy density (see, for example, Derrida jB*, Carmona and Hu 131). 
Recent work on a continuous model different to this can be found in Comets 
and Yoshida 0. 

In the next section we recall the framework which was developed in Jl] 
to extend some standard constructions from queueing theory to the context 
of geometric functionals of Brownian motion, and explain how this leads to 
the conjectured formula for the free energy density. Section |31 is devoted to 
the proof of the main result. In section |3] we prove that / is analytic and 
strictly convex, and record a large deviation principle that we will use in 
section |S1 to prove 

Acknowledgement: Research supported in part by Science Foundation Ire- 
land, Grant Number SFI 04/RP1/I512. The authors would also like to 
thank the anonymous referee for helpful suggestions which have led to a 
much improved presentation. 

2 Generalised Brownian queues 

In this section we recall the framework which was developed in Jl] to ex- 
tend some standard constructions from queueing theory to the context of 
geometric functionals of Brownian motion, and explain how this leads to the 
conjectured formula for the free energy density. 

The generalised Brownian queue is characterised as follows. Let B and 
C be two independent standard Brownian motions indexed by the entire 
real line, and write 



^{s,t) — Bt- Bs, C(s^t) — Ct- Cs 
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Fix m > and, for i G M, set 

r{t) = log / exp{S(s^t) + C^s,t) - m{t - s)} 

J —oo 

fis,t) = B^,^t)+r{s)-r{t) 
gis,t) = C^s,t) + r{s) - r{t) 

and define / : M ^ R by f{t) = f{0,t). 

To put this in context, the 'Brownian queue' is defined similarly but 
with 'log / exp' replaced by 'sup' (see, for example, Jl]). The usual M/M/1 
queue is defined similarly to the Brownian queue but with Brownian motions 
replaced by Poisson counting processes. Thus the Brownian motions Bt and 
mt — Ct can be thought of, respectively, as the arrivals and service processes, 
r as the queue-length process and / as the output^ or departure, process. 

In ^1] it is shown, using results of Matsumoto and Yor jl2l I13j . that 
the generalised Brownian queue is quasi-reversible, that is: / is a standard 
Brownian motion, and {/(s),s < t} is independent of r{t). We can thus 
consider a sequence of generalised Brownian queues in tandem and expect 
this 'queueing network' to have nice properties (analogous to the 'product- 
form solutions' of classical queueing theory). 

Let B,B^^\b'-^\... be a sequence of independent standard Brownian 
motions, each indexed by M, and let m > be a fixed constant. For — oo < 
s <t < oo, set 

/t 
ds exp{fi(^^t) + b'^1^^^^ - m{t - s)} 

his,t) = 5(,,t)+ri(s)-ri(t) 
g{s,t) = Sf^')^)+ri(s)-ri(t), 

for each A; = 2, 3, . . . set 

/t 
ds exp{/fc_i(s, t) + ^l^^^j^ - m{t - s)} 
-oo ' 

fk{s,t) = /fc_i(s,t) +rfc(s) -rfc(t), 

and for all k define //, : R ^ M by fk{t) = fk{0, t). 

Note that ri(t) is clearly stationary in t; to see that ri(0) < oo almost 
surely simply note that, with probability one, -B(s q) + -^(s o) ~'~ ^ ms/2 
for all s sufficiently negative (by Strassen's law of the iterated logarithm, 
for example). In fact, ri(0) has the same law as — logZ^, where Zm is 
gamma-distributed with parameter m: this is Dufresne's identity [Hl[7j. 

We first state the quasi-reversibility property, as presented in |14j . 

Theorem 1 1- fi and g are independent standard Brownian motions 
indexed by M. 
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2. For each t G M, g(s)), — oo < s <t} is independent o/{ri(s), s > 

t}. 

It follows from Theoremnthat ^1(0), r2(0), ... is a sequence of i.i.d. random 
variables, each distributed as — log Zm- By construction, we have 



" r ;>0 

^rfc(0)=log / d-u exp (5(^,^0) + 

1 L^— 00 



k=l 

X 



mu) 



^U<S1 <...<St,_1 <0 



(5) 



Applying the strong law of large numbers, it can be deduced (see |14| 
for details) that: 

Theorem 2 For each m > 0: 
1 f 

lim — log / dudsi . . . dsn-i 

n^OC n J -OD<U<Sl<...<Sn-l<0 

exp{mu + bI^^ , + ... + b')"'^ } = 

{u,si) (s„_i,0)J ^ ^ 

almost surely, where 

^(m) = ElogZrn = r'(m)/r(m) 
is the digamma function (and T is the Gamma function). 

Theorem 121 can be interpreted as follows. Let B denote the cr- field generated 
by the Brownian motions B^^\ B^'^\ . . . , and let ti, r2, ... be the points of 
a unit-rate Poisson process on IR+, independent of B. For t^^ti, . . . ,tn 
define 

En{to, ti,...,tn) = B^l^^^ + ... + ^(rj_^,j„) 

iC(io,ii,...,in-m+i) = expiB^;';. . + ... + B ' 



(iOjtl) ' ' ' itn-m,tn-m + l)' 



/? — 7?i 



By Brownian scaling, Theorem |21 is equivalent to: 

Theorem 3 For 9 ^ 0, 
1 



lim - logE[expi9En{0,n, . . . ,Tn))\B] = -2log\e\ - ^{1/9'), 

n^oo n 



almost surely. 
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Thus, if we set 



^(^^^ r -2iog|.l-^(i/n MO 

[ U, = 0, 

we have 

lim -logE[exp{eEniO,n, . . . ,Tn))\B] = Aie), 

n— >oo n 

almost surely. From the asymptotic expansion 

(6) 

k=l 

as X ^ OO (see, for example, Abramowitz and Stegun [I]), we have that A 
is finite and differentiable everywhere, with A(0) = A'(0) = 0. It follows 
that the sequence ^S„,(0, ri, . . . , r„) satisfies the following conditional large 
deviation principle: 

Theorem 4 Given B, ^En{0, ri, . . . , r„) satisfies a large deviation principle 
with good rate function 

A*(x) = sup[x9 - A{9)] 



almost surely. 

This is a quenched large deviation principle, associated with the condi- 
tional law of large numbers. For example. Theorem |3] implies that given B, 
y^En{0,Ti, . . . ,Tn) almost surely. Another implication is that for any 
X > 0, 

lim -logP{En{0,Ti, . . . ,Tn) > xn\B) = -A%x) 

n— ►oo n 

almost surely. For two other related large deviation principles see |14j . 

We will now describe how this relates to the Brownian directed poly- 
mer model. It is shown in Lemma IHl using Kingman's subadditive ergodic 
theorem that there exists a function 7 : R — > M such that given x < 0, 

lim i log / F„(xn,si, . . . ,s„_i,0)(isi . . = 7(x) (7) 

n^oo n Ja;n<si<...<s„_i<0 



almost surely, and it is shown in Lemma IHl that 7 is a concave function on 
(—00,0). Therefore, by Q, Theorem |21 and Laplace's method, we would 
expect 

—^{m) = sup[mx + 7(x)] = {—j)*{m) 



x<0 
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and hence by inversion, 7 = — (— The free energy density for our model 
of a directed polymer in a Brownian environment in 1 + 1 dimensions, defined 
in (j2j, can then be expressed in terms of the digamma function by first using 
the Brownian scaling property: 

hm -logZM = 7(-/?')-21og|/3| 

n—>oo n 

= -(-$)*(-/32)_21og|/3|. 

The heuristic argument above is made rigorous by the following Theo- 
rem, which is the main result of this paper. 

Theorem 5 Almost surely, 

hm i log = /(/?) 

where f is defined by Q, ^(m) = T'(m)/T(m) is the restriction of the 
digamma function to (0,00), and (— ^')* is the convex dual of the function 



3 Proofs 

We begin by defining for t > and x < 

Ln(t) = sup En{0,Si, . . . ,Sn-l,t) (8) 

0<si<...<s„_i<t 

K,A^) = log/ _ „,™)..,....»„ 

Zn{x) = / Fn{xn,si, . . . ,Sn-l,0)dSi . . .dSn-l (9) 

J Xn<Sl<...<Sn-l<0 

7„(x) = ^logZnix). (10) 
and recording the following lemma: 

Lemma 6 1. There exists a constant c G M such that ^Ln{nt) c\Jt 
almost surely and in expectation. 

2. There exists a function 7 : (— oo,0) M such that given x < 0, 

kQ^n{x)/n 7(x) (11) 

almost surely. 

3. The function 7 is continuous on (— cxd,0). 

4. Given x < 0, lim^^oo 7^(2;) = 7(x) almost surely. 
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Proof. The first part of the lemma fohows from Brownian scahng and 
Kingman's subadditive ergodic theorem (see for details) . For the second 
part, observe that 



ko,mix) + km,nix) = log/ Fn+l{xn, Si, . . . , Sn,0)dsi . . . dSn 

J xn<si<...<s„<0 

^n — m ^^'^^^n — m + 1 

and k is therefore super additive for fixed x. By construction we have the 
required conditions for Kingman's subadditive ergodic theorem, and so we 
may define a function 7 : (—00, 0) — > M by For the third part we have 
by Brownian scaling, for x < and 5 > x 

kQn{x-5) =d log/ exp (V-x + 5£'n+l(0, si, . . . ,Sn,n))dsi . . . dsn 

JO<si<... 

■ ..<Sn<n 

+n log(— X + 5) 

< log/ eX^p(^/^En+l{0,Sl,...,Sn,n))dSl...dSn 
JO<si<... 



.<Sn<n 



+{V-x + 6 - \/^)L„+i(n) + nlog(-x + 6) 
=d kQ,n{x) + + 5 - y/^)Ln+i{n) 

— nlog(— x) + nlog(— X + (5) (12) 

where =d denotes equality in distribution. Let c be the constant in the first 
part of this lemma; then by Brownian scaling and Slutsky's Lemma (see for 
example 0), ^L„+i(n) ^ c. By the second part of this lemma and (fT2|l . 

7(x — (5) < 7(x) + c(\/— X + 5 — a/— x) — log(— x) + log(— x + 5) 
By symmetry, 

K := inf £'„(0,si,...,Sn-i,n) =d -Ln(n) (13) 

and so we have similarly 



7(x — (5) > 7(x) — c(V— X + 5 — V—x) — log(— x) + log(— x + 6). 

So 

|7(x — 5) — 7(x)| < c\V—x + 5 — V—x\ 

+ |log(-x + (^) -log(-x)|. (14) 

For the fourth part we have 

fc0,n(x) < log ( e2^("'0'") / Fn{xn,S2,...,Sn,0)dsi...dSn) 

\ J Xn<Sl<...<Sn<0 / 

< log ( xne2^("'°'^) / Fn{xn,S2,...,Sn,0)ds2...ds 

\ J Xn<S2<...<Sn<0 



log (xne2^("'°'^)z„C 



X 
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where for y < x < 0, 

B{n,x,y)= max sup \B^'^\r,s)\. (15) 

i=n—l,n,n-\-l yn<r<is<xn 

From Borell's inequality and the Borel-Cantelh Lemma, there exists a null 
set M such that on its complement M'^, B(n, x, — > for all x, y, so 

7(x) < liminf 7„(x) a.s. 

n^oo 

Now for e > 0, 

koA^ + > log e"2^('"'°'''+') {x+e)n<si<xn Fn{xn, S2, ■ ■ ■ , Sn, 0)dsi . . . dSn 
\ J xn<S2<--- 

\ ...<Sn<0 

= logfene-2^("'°'^'+^)Z„(x) 



so 

7(x + e) > limsup7„(x) a.s. 



n— >oo 



and the result follows by the third part of this lemma. □ 
Let Q+ = Q n (0, oo), Q_ = -Q+. 

Lemma 7 There exists a null set M such that the following statement holds 
on its complement M^: 

lim„^oo7n(a::) = lix) and liminf„_^oo 7n(y) > 7{x) for every x G Q_ 
and y < X. 

Proof. Let x £ Q_ and y < x. From Q we have 

0)dsi . . . dSn-1 

J xn<si<...<s„-i<0 

= e (16) 
By Lemma El there exists a null set such that on M^, 

lim 7„(x) = 7(x). 

n— ►oo 

With J\f as in the proof of Lemma IHl let M = Af U Uxe(Q_ ■^x- ^ 
Lemma 8 Almost surely, lim„^oo "yn{x) = 7(x) for all x < 0. 
Proof. Choose x,y € Q_ with y < x. Then if z G [y, x], 

0)dsi . . . dSn-l 

J Zn<Sl<...<Sn-l<0 

= Ii{n, z) + hin, z) 
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where 

Ii{n,z) = e'°(^"'^")Z„(x) 



h{n,z) = Fn-iisi,...,Sn-i,0)ds2...dsn-ie (^"'»i)(isi 



Now Ii{n,z) < e^("'^'2')Z„(x) and 

Vn — 1/ 

therefore 



sup Zn{z) < e^("'"'^)^„(x) + (x - y)ne2^("'^'J')Z„_i(^). (17) 



2G[j/,a;] 1 

Take as in the proof of Lemma [7| Now 

Vn — 1/ 



so hminf„^oo ^log2:„_i(^^j > -f{y) on M""; and if e G then 



j/n<S2<---<s„_i<0 



^^-2B{n,0,y+e) 



— 1 
hence 

7(y + e) > hmsup-logZ„_i - 

n Vn — 1 



n^oo 



on A^"^, and letting e ^ in Q gives hm^^oo ^logZ„_i(^^j = 7(y) on 
A^*^. Then by 1)17^ and the proof of Lemma [3 on A^'^ 

1 1 
7(2;) < hminf — log inf ^^^(z) < limsup — log sup Zn{z) 

< limsup - log e^("'^'S')Z„(x) V 

n^oo n 

limsup - log(x - y)ne2^("'^'^)Z„_i (-^ 
= 7(2;) V7(y) =7(?/)- (18) 
The result now follows from 1)141) . □ 



9 



Lemma 9 The function 7 is concave. 
Proof. For x, y < and a G (0, 1), 



7n(ay + (1 - a)x) > ^^Gn + -H^ (19) 

n n 



where 



= il°g/.y+(l-a).)n<.,<...^n"((«y + (l-«)^K«l, 
•■■<S[c«n]-l<(l-o)a;n 

• • • , S[an]-1, (1 - a)a;n)(isi . . . dS[an]-l 

kn = n — [an] + 1 

loe [, , Fu((l — a)xn, s\ry„], . . . 

...<s„_i<0 

. . . , 0)(iS[an] • • • dSn-l 



Hn 



Now 



/ an \ 

/(l-a)n , 
-?^n = 7fc„ ^ a; . 



Choose w,u £ Q_ with w < x < u, and choose e, 5 > 0. Then 3 no G N 
such that y n > no, w < ^ < ^- Also, by (fTH|) 3 ni € N such that 

V n > ni, 

P(j{u) - e < - log inf Zn{z) < - log sup Z„(z) < 7(t(;) + e) > 1 - J 
V n 2e[ii),u] n ^ 

Therefore if A:„ > no V ni, 

P(7(n) -e<Hn< i{w) + e) > I - 5 

and therefore ff„ =^ 7(2;) • Similarly G„ =^ 7(y), and hence by Slutsky's 
Theorem 

+ aj{y) + (1 - a)7(x) 

n n 

Hence by (fT9|) . 

7(0;^ + (1 - q;)x) > Q;7(y) + (1 - Q;)7(a;). 

□ 
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Proof of Theorem Choose m > 0. Define a probability density 
function on (— cxd,0) by 

K^\dx) = ^— — r- exp n(mx + 7„(x))(ix 

where 



H„(m) = / exp n(m2; + 7„(2;))(ix. 

J — oo 

In the nomenclature of statistical physics, icj^"^ is the Kac density and 
H„(m) is the grandcanonical partition function. Theorem |^ says that given 
6 > —m, the convergences 

-logH„(m) -^(m) (20) 

n 

1 1-0 

- log / e"^^IK(™)(da;) ^ ^(m) - ^(m + 0) =: A^f^) 
n J-oo 

hold almost surely as n ^ cxd. Choose e > 0; then since AJ^(A^(0) + e) and 
A5^(A^(0) — e) are strictly positive, we can apply the Chernoff bound to give 
]K^i™'^((A^(0) — e, A^(0) + e)) ^ 1 almost surely as n ^ oo. Letting e ^ 
in Q gives that K^™''' is almost surely concentrated on A^(0) = — ^''(m) as 
n ^ oo. Therefore for any x G Q_ and e > 0, using (|16j) we have 



H„(m) > / expn(my + 7„(y))dy 

J x—e 

px 

> exp{n{m{x - e) + 7„(x))} / exp B^^^l^^^^dy 

J x~e ' 

> e exp{n(m(x — e) + ^nix)) — B{n,x,x — e)} 
Therefore by Lemma El and (|2()jl . 

— ^'(m) > ?7i(x — e) + 7(x) 

and we may let e — > and appeal to the regularity of ^ to conclude that for 
all X < 

7(x) < inf (-mx - *(m)) = -(-^)*(x). 

mGQ+ 

For the reverse inequality we note that for x G Q_ and e G (0, — x) we may 
choose m > such that — ^''(m) = x + e (see for example m2])- Then by 
(CSl), 

]K(rH(x,x + 2e)) = r^'%xp{n(my + 7„(y))}dy 

< 7— — - exp{n(m(x + 2e) + 7„(x)) - B(n, x + 2e, x)} 
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Therefore using Lemma |H1 

= lim -log]K(7')((x,x + 2e)) < m(x + 2e) +7(x) + ^'(m) 

n— >oo n 

almost surely, in which case we may let e ^ to conclude that 
7(x) > inf [—mx — ^{m)) = —{—'^)*{x). 

Therefore 7(x) = — (— Finally the last part of Lemma IHl and 
Brownian scaling allow us to conclude that almost surely, 

n-^oo n \^ i- : p = U 

as required. □ 

4 Analyticity of the free energy density and a large 
deviation principle 

Theorem 10 The function f defined in ^ is analytic and strictly convex 
on R, /'(O) = 0, and lim^^^ /(/?)//? = 2. 

Proof For x < we have 

= sup[xe + "^{9)] 

Denoting by the nth derivative of the function and noting that ^2 is 
strictly negative everywhere (see for example [2j), we have 

and since is an invertible analytic function with nonzero derivative, its 
inverse is analytic. Therefore / is analytic everywhere except possibly at 0. 
To investigate the behaviour of / near 0, let a = ^tH/S^)- Then 

f{(3) = a^i{a)-^{a)-log^i{a) (21) 

Now a — > 00 as /? — > 0, and from we have 

+ ^ + ^ (22) 

therefore recalling 

/(/?) = l + 0(a-i) (13^0) 
One further application of H22() to this expression gives /'(O) = 0. 



12 



Since / — 1 is an asymptotic logarithmic moment generating function 
(see the proof of Lemma^Jbelow), / is convex. The regularity of ^ implies 
further that / is strictly convex. 

The power series 

oo 

logr(l + z) = - log(l + z) + z{l -0 + JZi-mCin) - l]z^/n 

n=2 

where ^ is Euler's constant and Q is the Riemann Zeta Function, is valid for 
\z\ < 2 therefore, since '^{z) = ^logr(z), it may be differentiated to 
give power series for ^ and ^'i in a neighbourhood of the origin. Substituting 
in (|2T|) and letting a — > 0, we obtain lim^j^oo fiP)/P = 2. □ 

Lemma 11 Let ti < . . . < t„_i be the order statistics for n — 1 indepen- 
dent random variables having the uniform distribution on the interval [0, n] . 
Almost surely, conditional on B, the random variable ^£J„(0, ri, . . . , rn_i, n) 
satisfies a large deviation principle with rate function (/ — 1)*. 

Proof. Choose /? G M; then by (UJ 

E[exp(/3K(0,ri,...,r„_i,n))|^] = i-^Z„(/3) 
hence by Stirling's formula and Theorem |S1 

hm -log^[exp(/3^„(0,ri,...,r„_i,n))|e] = /(/5)-l (23) 

n— ►oo n 

almost surely. Now \i a < u then Zn{a) < e^^^~"''>'^" Zn{v), where Vn was 
defined in ^T^: hence for a, /3 G Q with (3 > a 

f{a) < liminf — log inf Z„(z^) + (/3 — a)c 

n-^oo n ue{a,/3) 

< limsup-log sup Z„(i^) + (/3 - a)c </(/?) + 2(/3 - a)c 

almost surely, where c was defined in Lemma IHl Therefore by the continuity 
of /, there exists a null set J\f such that on M'^, the convergence in (|23|) holds 
for ah /3 G M. □ 

5 Connection with random matrices and a Brow- 
nian directed percolation problem 

It was shown in jllj that with Ln defined as in Q, for each t >0, 

lim —Ln{nt) = 2\ft 

n—>oo n 
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almost surely. In the notation of Lemma El this says that c = 2. This 
result can also be deduced from random matrix theory using the fact 
that Ln{l) has the same law as the largest eigenvalue of a n x n GUE 
random matrix. We note here that the present work allows us to deduce the 
inequality c > 2. Since Ln{n) > En{0, ti, . . . , tn-i,n) we have that given /3, 

-Ln{n) > ^-log£;[exp(/3K(0,ri,...,T„_i,n))|^] 

n pn 

almost surely. Letting n ^ oo and then /3 — > oo, using ()23|) and Theorem 
IIUI gives the required inequality. 
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